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Abstract. The connection between helically isotropic MHD turbulence and mean-field 
dynamo theory is reviewed. The nonlinearity in the mean-field theory is not yet well 
established, but detailed comparison with simulations begin to help select viable forms 
of the nonlinearity. The crucial discriminant is the magnetic helicity, which is known 
to evolve only on a slow resistive time scale in the limit of large magnetic Reynolds 
number. Particular emphasis is put on the possibility of memory effects, which means 
that an additional explicitly time-dependent equation for the nonlinearity is solved 
simultaneously with the mean-field equations. This approach leads to better agreement 
with the simulations, while it would also produce more favorable agreement between 
models and stellar dynamos. 



1 Introduction 

In an early paper Parker [0 identified cyclonic convection as a key process for 
converting large scale toroidal magnetic field into poloidal fields that have co- 
herence over about half a hemisphere. This process is now generally referred to 
as the a-effect, although it may arise not only from thermal buoyancy |2 , but 
also from magnetic buoyancy fi P]5|] , the magneto- rotational instability ||,7|], or 
some other magnetic instability |J. In each case the effect of the Coriolis force 
together with some kind of radial stratification is crucial for making the motions 
helical Upward moving fluid expands, and the Coriolis force makes it rotate 
retrograde, causing negative (positive) kinetic helicity in the northern (southern) 
hemisphere. Downward moving fluid contracts, rotates in the prograde direction 
and contributes in the same sense to negative (positive) kinetic helicity on the 
northern (southern) hemisphere. This causes a positive a-effect in the northern 
hemisphere, but if magnetic stresses and shear become strong (for example in 
accretion discs) the sign may reverse [ p~0[JTT[ ] . 

When combined with differential rotation, the main outcome of a-effect mod- 
els is the possibility of cyclic magnetic fields associated with latitudinal migra- 
tion. The first global (two-dimensional) models were presented by Steenbeck & 
Krause ]T^|, but similar models, with different physics, are still being studied 
today |13| , [14|Jl5|| . The migratory behavior is best seen in contours of the longi- 
tudinally averaged mean magnetic field versus latitude and time, which should 
show tilted structures converging to the equator. Such plots can be compared 
with the solar butterfly diagram of sunspot numbers (so called because the struc- 
tures resemble a sequence of butterflies). 
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A key property of all these models is that not only the motions are helical, 
but the large scale magnetic field itself is also helical. Of course, not all dynamos 
require helicity, but nonhelical dynamos tend to generate preferentially small- 
scale fields p6[ . In a recent attempt, Vishniac & Cho |17 proposed a mechanism 
relevant in particular to accretion discs where shear is strong. Their mechanism 
would not lead to the production of net magnetic helicity, but numerical simula- 
tions |l8| failed so far to show convincingly large scale dynamo action based on 
the proposed mechanism. Shear does produce large scale fields, but only in the 
toroidal direction. It does not explain the latitudinal coherence of the field over 
several tens of degrees (corresponding to several hundred megameters). On the 
other hand, there is direct observational evidence that the solar magnetic field 
is indeed helical. (We shall return to observations in Sect. ||.) 

The trouble with helical fields is that magnetic helicity is conserved by the 
induction equation in the ideal limit and can only change on a resistive time 
scale, provided there is no significant loss through boundaries (at the surface 
or the equator, for example). This approximate magnetic helicity conservation 
leads to magnetic field saturation on a resistive time scale [Q. Depending on 
how effectively the boundaries transmit magnetic energy and helicity, the final 
saturation amplitude will be lowered if losses occur preferentially on large scales 
while the (linear) growth rate of magnetic energy (past initial saturation) re- 
mains unchanged p(| . In this sense final saturation can be achieved earlier. The 
above results are particularly clear when the flow is nearly fully helical, i.e. when 
the normalized kinetic helicity, ef = (<*> • u)/(w rms it rms ), where u> = V x u is the 
vorticity, is large. In the sun, and probably in all other celestial bodies with 
rotating turbulence, the relative kinetic helicity is small; ef ~ 5%. It is however 
this small helical fraction of the turbulence that is responsible for the a finite 
but small a-effect, so a proper understanding of its dynamics is crucial if one 
wants to build models based on the a-effect. Below we shall also discuss the 
alternative that astrophysical dynamos may shed preferentially small scale he- 
lical fields through the boundaries. This could theoretically enhance large scale 
dynamo action pil , |2^ , p3| . 



2 Magnetic helicity production 

Before we begin discussing the magnetic helicity problem and possible remedies 
we need to be sure that the solar dynamo is indeed likely to involve significant 
amounts of magnetic helicity. There is direct observational evidence that the field 
of the sun is actually helical. Firstly, active regions are known to have systemati- 
cally different signs of current helicity in the two hemispheres p^ , p5p^ ,p7| ; pref- 
erentially negative (positive) in the northern (southern) hemisphere. Secondly, 
magnetic helicity flux from the solar surface has also been inferred and this con- 
firms the same sign as that of the current helicity. The magnetic helicity flux 
driven by the surface differential rotation has been estimated by Berger & Ruz- 
maikin pq ], who find a total magnetic helicity flux on the order of 4 x 10 46 Mx 2 
over the 22 year solar cycle. Similar values were also found by DeVore p9|. 
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Finally, Chae |30| estimated the magnetic helicity flux based on counting the 
crossings of pairs of flux tubes. Combined with the assumption that two nearly 
aligned flux tubes are nearly parallel (rather than antiparallcl) his results again 
suggest that the magnetic helicity is preferentially negative (positive) in the 
northern (southern) hemisphere. 

The magnetic helicity is noisy, i.e. the sign can fluctuate and has only on 
average systematic behavior. This reflects the fact that only a fraction of the 
turbulence is helical. Thus, detailed understanding and measurements of the 
departures from systematic behavior is just as important as understanding and 
recording the systematic behavior. 

In astrophysical flows, kinetic helicity can be generated in rotating strati- 
fied turbulence ||. Such flows are intrinsically anisotropic. Whilst this is not a 
problem for numerical simulations, it definitely complicates the theoretical un- 
derstanding and one should not be surprised if some fundamental aspects of 
mean-field theory (e.g. a proportional to — (u> • u}) are not recovered. We just 
mention that under certain conditions, stratified rotating flows exhibit an ex- 
tensor whose vertical (z-z) component has the opposite sign as the horizontal 
(x-x and y-y) components [ |32|]3^ , |34| , |35| |.) 

In the following we concentrate on the isotropic helical aspects of the tur- 
bulence. This is accomplished by driving the flow with random polarized waves 
in a periodic domain. In most of the cases we use fully helical turbulence, but 
in many estimates the fraction of helicity enters only as an additional scaling 
factor Q . The main goal here is a better understanding of the a-effect and the 
nonlinear feedback when the field becomes dynamically important. We there- 
fore discuss in detail the perhaps simplest possible system: the a 2 -dynamo in a 
periodic box. 

In spherical geometry, the term a 2 -dynamo refers to the fact that both large 
scale poloidal and toroidal fields are maintained against ohmic decay by the a- 
effect. By contrast, the a/2-dynamo is one where the large scale toroidal field is 
generated mostly by differential rotation (the f2-effect) and the a-effect can be 
neglected by comparison. If the a-effect is not neglected one speaks of an a 2 Q- 
dynamo. We stress that the a 2 -dynamo has nothing to do with the so-called 
small scale dynamo. These are turbulent dynamos operating only on scales less 
than the energy-carrying scale of the turbulence. They are quite common if the 
flows arc non-helical. By contrast, both a 2 and aJ7-dynamos also generate fields 
on large scales, but they are necessarily accompanied by some level of small scale 
fields as well. 

In its simplest form the a 2 -dynamo equations for isotropic a and turbulent 
diffusivity rjt can be written as 



— = V x (aB - 77TA*oJ) , (1) 

where J = V x B//xo is the mean current density and /iq is the magnetic 
permeability. This equation permits plane wave solutions of the form B = 
B exp(At + ifc • x), with the dispersion relation A = ±|a|fe — ri^k 2 , where k — \k\. 
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The maximum of A is where dA/dfc = 0, which yields 

k = U = a/(2rn). (2) 

For a periodic box of size L 3 , the most easily excited mode has k — k\ = 2w/L, 
and the k vector can point in any of the three coordinate directions. For a < 
(the case considered in Ref. the three possible eigenfunctions are 

(cosfcizA / \ / sinfeiyN 

sinfciz , B W =B cosfcii , B (y) = B \ , (3) 
I \smkix J \coskxy J 

where we have ignored arbitrary phase shifts in any of the three directions. All 
three solutions have been found in the simulations p9| . 

In the simulations there is of course no explicit a-cffcct in the usual sense, be- 
cause we just solve the primitive MHD equations. The turbulence does, however, 
display collective behavior - just as it is expected based on mean-field a 2 -dynamo 
theory, as explained above. We begin with the simulations. 



3 Helical turbulence: prototype of an CK 2 -dynamo 

We consider a compressible isothermal gas with constant sound speed c s , con- 
stant dynamical viscosity p, and constant magnetic diffusivity r/. To make sure 
the magnetic field stays solenoidal, i.e. V • B = 0, we express B in terms of the 
magnetic vector potential A, so the field is written as B = V x A. The governing 
equations for density p, velocity u, and magnetic vector potential A, are given 

by 

D In p _ , , . 

-DT = -V-u, (4) 

^ = -c 2 V hxp + — + ^(V 2 u + i VV • u) + /, (5) 

Dt p p 3 

dA 

— = u x B - rip J - V<f>, (6) 

where D/Dt = d/dt + u ■ V is the advective derivative. The current density, 
J = V x B//i , is obtained in the form p 3 = -V 2 A + VV • A. We often use 
cj) = as a convenient gauge for the electrostatic potential. Other frequent choices 
are 4> = —rfV ■ A, <f> = u- A, or combinations of these j36|. The Coulomb gauge, 
V-A = 0, corresponds to (j> = — V-E, where E = — uxH+rjp J, but the original 
reason for solving for A instead of B was just to get rid of the solenoidality 
condition, so the Coulomb condition has computational disadvantages. 

For the following it is useful to recall that each vector field can be decom- 
posed into a solenoidal and two vortical parts with positive and negative helicity, 
respectively. These are also referred to as Chandrasekhar-Kendall functions. Al- 
though it is often useful to decompose the magnetic field into positive and nega- 
tive helical parts, we also use the helical fields (with positive helicity) as forcing 
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function / of the flow. We restrict ourselves to functions selected from a finite 
band of wavenumbers around the wavenumber kf , but direction and amplitude 
are chosen randomly at each timestep. Details can be found in Ref. ]l9| . Similar 
work was first carried out by Meneguzzi et al. p7j , but at the time one was barely 
able to run even until saturation. Throughout the nineties, work has been done 
on forced ABC flows (3^j3^,|{J. In none of these investigations, however, the is- 
sue of resistively slow magnetic helicity evolution past initial saturation has been 
noted. It is exactly this aspect that has now become so crucial in understanding 
the saturation behavior of nonlinear dynamos. We begin by discussing first the 
linear (kinematic) evolution of the magnetic field. 

3.1 Linear behavior 

If the magnetic Reynolds number, defined here as R m — u ms /(r]k{), exceeds a 
certain critical value, Rm Tlt \ there is dynamo action. For helical flows, the so 
defined -Rm"*' 1 is between 1 and 2 (see Table 1 of Ref. jl9]]). In the supercritical 
case, 

Rm > Rm it] , the field grows exponentially with growth rate A, which 
scales with the inverse turnover time, u rms /k{. The resistively limited saturation 
behavior that will be discussed below in full detail has no obvious correspondence 
in the kinematic stage when the field is too weak to affect the motions by the 
Lorentz force jyj. Nevertheless, there is actually a subtle effect on the shape 
of the eigenfunction as R m increases. Before we can appreciate this, we need to 
discuss the effect the kinetic helicity has on the field. 

A helical velocity tends to drive helicity in the magnetic field as well, but 
in the nonresistive limit magnetic helicity conservation dictates that (A • B) = 
const = if the initial field (or at least its helicity) was infinitesimally weak. 
Thus, there must be some kind of magnetic helicity cancelation. Under homoge- 
neous isotropic conditions there cannot be a spatial segregation in positive and 
negative helical parts. Instead, there is a spectral segregation: there is a bump 
at the forcing wavenumber and another 'secondary' bump at somewhat smaller 
wavenumber. The two bumps have opposite sign of magnetic helicity such that 
the net magnetic helicity is zero. At the forcing wavenumber, the sign of mag- 
netic helicity agrees with that of the kinetic helicity, but at smaller wavenumber 
the sign of magnetic helicity is opposite. At small R m , this secondary peak can be 
identified with the wavenumber where the corresponding a 2 -dynamo has max- 
imum growth rate; see Eq. (§). Simulations show that as R m increases, fc max 
approaches ^kf [p3| . This agrees qualitatively with earlier results |4l];^2| which 
suggested that the magnetic helicity approaches zero in the high-i? m limit. 

3.2 Nonlinear behavior 

Eventually, the magnetic energy stops increasing exponentially. This is due to 
the nonlinear terms, in particular the Lorentz force J x B in Eq. (jHJ), which begins 
to affect the velocity field. The temporal growth of the power spectra saturates, 
but only partially; see Fig. [I], where we show data from a run with forcing at 
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Fig. 1. Power spectra of magnetic energy of Run 3 of Ref. |19[. During the initial 
growth phase the field saturates at small scales first and only later at large scales (left 
hand panel). Later, when also the large scale field saturates, the field at intermediate 
scales (k = 2, 3, and 4) becomes suppressed. In the second panel, intermediate times 
are shown as dotted lines, t = 700 is shown in solid and t = 1600 is shown as a thick 
solid line. The forcing wavenumber is kt = 5. 



wavenumber fcf — 5. In the left hand panel we see that by the time t = 600 
the power spectra have saturated at larger wavenumbers, k <; 3. It takes until 
t ~ 1600 for the power spectra to be saturated also at k = 1 (right hand panel of 
Fig. [j]). In order to see more clearly the behavior at large scales, we show in Fig. || 
data from a run with fcf — 27 and compare spectra in the linear and nonlinear 
regimes. In the linear regime, all spectra are just shifted along the ordinate, so 
the spectra have been compensated by the factor Mj n ; exp(2Ai), where A is the 
growth rate and -Mi n ; the initial magnetic energy. In the nonlinear regime the 
bump on the right stays at approximately the same wavenumber (the forcing 
wavenumber), while the bump on the left propagates gradually further to the 
left. As it does so, and since the amplitude of the secondary peak even increases 
slightly, the net magnetic helicity inevitably increases (or rather becomes more 
negative in the present case). But because of the asymptotic magnetic helicity 
conservation, this can only happen on a slow resistive time scale. This leads to 
the appearance of a (resistively) slow saturation phase past the initial saturation; 
see Fig. ||. 

3.3 The final saturation value 

In a periodic box with helically driven turbulence, the final saturation value of 
the magnetic field is determined by the ratio of the size of the domain to the scale 
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Fig. 2. Power spectra of magnetic energy of positively and negatively polarized parts 
(M^ and Mr) in the linear and nonlinear regimes. The spectra in the linear regime 
have been compensated by the exponential growth factor to make them collapse on 
top of each other. Here the forcing wavenumber is in the dissipative subrange, kf = 27, 
but this allows enough scale separation to see the inverse transfer of magnetic energy 
to smaller k. The data are from Run B of Ref. vAM. 



of the forcing. This is best seen by considering the magnetic helicity equation 

A( A -B)=-2 W o(J-B). (7) 

If (A ■ B) were not gauge invariant (for example if there are open boundaries), 
Eq. (0) would be useless. In particular, (A • B) will in general not be constant in 
the steady state (see Fig. 2 of Ref. (^3), for an example). One therefore has to go 
to the gauge-invariant relative magnetic helicity of Berger & Field . This has 
been done in Refs ^,^3|. In the present case of periodic boundaries, however, 
(A • B) is automatically gauge invariant and therefore a physically meaningful 
quantity, so it must be constant in the steady state. Equation (|?]) says that then 
the current helicity, (J • B), must vanish. At first glance this seems to be in 
conflict with the idea of building up a helical large scale field. The solution is 
that there must then be an equal amount of small scale current helicity so that 

(J -B) = (J -B) + (j -b) = (in the steady state). (8) 

2 

If the field is fully helical, we have /.to (J • B) = =Ffci(B ) and /io(j • b) = ±fcf (b 2 ), 
where the upper (lower) sign applies to the case where the small scale helicity 
at the forcing scale is positive (negative). We then have from Eq. (^) 

<B 2 ) = (fc f Ai)<b 2 ). (9) 
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Fig. 3. The three stages of the magnetic field growth: exponential growth until initial 
saturation (when (B 2 )//io = (pu 2 )), followed by a (resistively) slow saturation phase. 

2 

In this plot we have used po = 1. The energy of the large scale magnetic field, (B }, is 
shown for comparison. The data are from Run 3 of Ref. . 

To a good approximation, (b 2 ) 1 / 2 will be close to the equipartition field strength, 

-Bcq, so 

(B 2 )/S 2 q a kt/ki > 1. (10) 

This means that the energy of the large scale field must, in the final state, be 
in super-equipartition by a factor approximately equal to the scale separation. 
We recall that this applies to the case of periodic boundaries. For closed (e.g. 

2 

perfectly conducting) boundaries, (B } can be even larger than kf/k\ times the 
equipartition value pj|. This is because the large scale field is no longer fully 
helical, while the small scale field still is. In the presence of open boundaries, on 
the other hand, the large scale field will generally be smaller than suggested by 
Eq. (|Io|), unless the boundaries transmit preferentially small scale fields (Sect. ||). 

3.4 Sensitivity to using hyperdiffusivity 

The above statements can readily be generalized to the case where the usual 
magnetic diffusion operator, 77V 2 , is replaced by hyperdiffusion, (— l) n ~ 1 T] n W 2n 
with n = 2. This implies that the diffusion has now become more strongly 
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wavenumber dependent; from r\k 2 to 772A: 4 . If the diffusion is the same at small 
scales, then the diffusion at large scales must be significantly smaller in the 
hyperdiffusivc runs. This leads to a dramatic increase of the resistive saturation 
time. At the same time the final saturation field strength is no longer given by 
Eq. (|To|). The rate of magnetic helicity dissipation is now no longer proportional 
to k, but to fc 3 . Therefore the final saturation field strength is given by 

(B 2 )/S c 2 q a (Vfci) 3 » I- (11) 

This result was confirmed also numerically jf4|. Again, this applies to periodic 
boundaries. Hyperdiffusion has been used in the past in connection with open 
boundaries ||(l5], but it is not clear how serious the possible artifacts from 
hyperdiffusion would be in such cases with open boundaries. 



3.5 The magnetic helicity constraint 

The case of periodic boundaries is particularly useful as a benchmark to all 

dynamos exhibiting large scale dynamo action due to the helicity effect. Here we 

2 

discuss the functional form B (t) for the late saturation phase. 

Equation (^) allows us not only to determine the final saturation strength, 
but also the approximate time evolution near saturation. As before, we make 
the assumption of fully helical fields. However, given that prior to saturation 
I (j • b) I « I (J ■ B) I , we must have | (a ■ b) | <C | (A ■ B) | , so we can set 

(A-B)«(A.B)=Tfcr 1 (B 2 >, (12) 

where the upper (lower) sign refers to positive (negative) small scale (kinetic 
and magnetic) helicity. Inserting this into Eq. (Q) we have 

fcf 1 A(B 2 ) = -2 v h{B 2 ) + 2 v k f {b 2 }. (13) 

The small scale field saturates first, so Eq. ( |l3| ) can then be integrated to get 
the subsequent evolution of (B ) toward full saturation |T^] 



kt 



(B-) = -^<b*) i- e -^(t-u at ) (for i >tsat)] (14 ) 
k\ L 

where i sa t is the time when the small scale field has reached saturation. Equa- 
tion (]l4|) is what we usually mean by the magnetic helicity constraint. 



3.6 Inverse cascade versus a-effect 

The process outlined above can be interpreted in two different ways: inverse 
cascade of magnetic helicity and/or a-effect. The two are similar in that they 
tend to produce magnetic energy at scales larger than the energy-carrying scale 
of the turbulence. As can be seen from Figs [I] and @, the present simulations 
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support the notion of nonlocal inverse transfer |TS| . This is not really an inverse 
cascade in the usual sense, because there is no sustained flux of energy through 
wavenumber space as in the direct Kolmogorov cascade. Instead, there is just a 
bump traveling to smaller k in wavenumber space. In that respect, the present 
simulations seem to differ from the Eddy Damped Quasi-Normal Markovian 
(EDQNM) closure approximation p6j . 

The other interpretation is in terms of the a-effect. We recall that there is a 
wavenumber fc max where the growth of the large scale held is largest; see Eq. (||). 
For reasonable estimates, /c max coincides with the position of the secondary bump 
in the spectrum (Ref. Sect. 3.5). This can be taken as evidence in favor of 
the a-cffect. In the nonlinear regime, the secondary bump travels to the left 
in the spectrum (i.e. toward smaller k). In the EDQNM picture this has to do 
with the equilibration of kinetic and current helicities at progressively smaller 
wavenumbers, which then saturates further growth at that wavenumber, but per- 
mits further growth at smaller wavenumbers until equilibration occurs. Another 
interpretation is simply that if a is quenched to a smaller value, /c max = a/(2rfr) 
peaks at smaller wavenumbers where the growth has not yet saturated, until 
equilibration is attained also at that scale. 

3.7 Implications 

The growth of the large scale magnetic field can be interpreted as being due 
to the a-effect. Consequently, a slow-down in the final saturation phase must 
have to do with a suppression of a. According to closure models [Q, the a- 
effect is really the residual between two competing effects: a kinematic helicity 
effect (which itself decreases somewhat for strong magnetic fields), and a current 
helicity effect of opposite sign. In Sect. ^ we shall use this phenomenology in a 
mean-field model where an explicitly time-dependent equation for the current 
helicity is solved. It turns out that in this model the late saturation phase is 
resistively limited-just like in the simulations. 

It is somewhat worrisome that in the nonlinear regime the value of a depends 
on the microscopic magnetic diffusivity, which is very small in most astrophysical 
situations. One might therefore be concerned that in astrophysical dynamos the 
saturation would be unacceptably slow. In order to say more about a and also 
the turbulent magnetic diffusivity rj t , we need to determine how a and rj t depend 
on B. This can either be done directly jl7|, which is difficult and the results are 
noisy, or we can compare with models that incorporate the effects of a and 7ft 
quenching. Before we do this (Sects || and 0), we first want to assess the effects 
of boundaries. 

4 Open boundaries: good or bad? 

Boundaries generally lead to a loss of magnetic field both on small and large 
scales. Losses at large scale tend to lower the saturation field strength of the 
mean field. Losses of small scale fields can, at least in principle, enhance the 
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large scale field p4p2[ . This has been demonstrated in an idealized numerical 
experiment [ p3[ , where the magnetic field at the forcing wavenumber and beyond 
had been removed in regular time intervals; see Fig. ^. We discuss this now in 
more detail. 
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Fig. 4. The effect of removing small scale magnetic energy in regular time intervals At 
on the evolution of the large scale field (solid lines). The dashed line gives the evolution 
of (B ) for Run 3 of Ref. jl!^ (where no such energy removal was included) in units 
of _Bc q = popo{u 2 ). The two solid lines show the evolution of (B ) after restarting the 
simulation from Run 3 of Ref. jl9| at At = 20 and At = 80. Time is scaled with the 
kinematic growth rate A. The curves labeled (a) give the result for At = 0.12A -1 and 
those labeled (b) for At = 0.4A -1 . The second panel shows, for a short time interval, 
the sudden drop and subsequent recovery of the total (small and large scale) magnetic 
energy in regular time intervals. (Adapted from Ref. [E3j.) 



4.1 Enhancement through losses at small scales 

The enhancement of large scale field by losses at small scales may seem some- 
what mysterious. One way to interpret this result is by saying that the slow 
growth occurred because the energy of the small scale magnetic field has already 
reached the level of the kinetic energy, and only the large scale field has not 
yet saturated. After small scale magnetic fields have been removed, the field is 
for a short time interval in sub-equipartition at small scales and so the overall 
field (both at small and large scales) can then grow further during the short 
time interval during which the small scale field has not yet fully recovered to 
the equipartition value. The effect of a single such event is small, but many such 
events can produce a significant effect. This is exactly what is seen. Another way 
of interpreting this result is in terms of mean-field theory where the a-effect is 
saturated by a cancelation of kinetic and current helicities. If small scale mag- 
netic fields are removed, the residual a-effect can be larger for some time interval, 
which then allows the field to grow further. In the following we illuminate this 
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result further by considering a modified magnetic helicity constraint for the case 
of open boundaries. 



4.2 The modified magnetic helicity constraint 



In Sect. 3.5 we have discussed an equation for the evolution of the magnetic 
energy of the large scale field at late times. Here we have assumed that there is 
no loss of magnetic energy and magnetic helicity through the boundaries. This 
equation has been generalized to account for losses of large scale magnetic helicity 
[ p0|]23|] . The idea is that there will be a magnetic helicity flux that is proportional 
to the gradient of the large scale magnetic helicity density (in a fixed gauge) , and 
hence to the gradient of the magnetic energy density. This gives rise to an extra 
diffusion term, and hence to a renormalized, effective magnetic diffusivity, rj^j? , 

i.e. the term 2r]ki (B 2 ) has to be replaced by 2ry^g ki(B 2 ). Therefore, Eq. ( p"3| ) 
takes the form |E3| 



A ; r 1 ^(B 2 ) = -2^ ) fc 1 (B 2 ) 



2^ f (b 2 }, 



which has the solution 

-2, T)k{ 



<B*> 



(b 2 ) 



l 



-2T)$fef(t-t 8 »t) 



(for t > taat), 



(15) 



(16) 



Note that the saturation amplitude is decreased by a factor r)/r)\jx compared with 

Eq. (p^), but at the same time the e-folding time has decreased to (2»7^jj A^) -1 . 
We return to this behavior in the next subsection. 

When the losses through the surface occur preferentially at small scales, an 
effective diffusivity would instead affect the small scale helicity flux. Therefore, 
Eq. ( |l4| ) takes then the form 

d_ 

di' 

which has the solution 



K 1 ^) = -2 v h{B 2 ) + 2r ? «fc f (b 2 ), 



(17) 



(B 2 ) = 



%B K i /v2 



1 _ e -2i7*?(t-W) 



(for t > i sat ), 



(18) 



Note that the saturation amplitude is now increased by a factor ry c g jr\ compared 
with Eq. (|lj), but the e-folding time, {2rjk\)~ 1 , is still resistively limited. This is 
in good with what is seen in the simulations; see Fig. ||. In reality, there will be 
both small and large scale losses, so the large scale magnetic energy is expected 
to evolve according to 



(f)i, 



1 _ e -2»7S *!?(*- W) 



(for t > t^). 



(19) 



This equation allows time scales and saturation amplitudes that are not resis- 
tively limited. 
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4.3 Simulations with open boundaries 

So far, simulations have not yet shown that the losses of small scale magnetic 
fields are actually stronger than those of large scales fields. Simulations with 
a vertical field (pseudo-vacuum) boundary condition have shown that most of 
the magnetic energy is lost at small scales |2(|. The way this affects the slow 
resistively limited saturation process discussed earlier is by simply cutting off 
the saturation process at an earlier time, without changing the approximately 
linear slope past the initial saturation; cf. Eq. (|l6|). In Fig. || we demonstrate 
a very similar behavior in another system which is actually periodic, but the 
helicity of the forcing is modulated in the z-direction such that the sign of the 
kinetic helicity changes in the middle. One can therefore view this system as two 
subsystems with a boundary in between. This boundary would correspond to 
the equator in a star or the midplane in a disc. 




Fig. 5. Evolution of the magnetic energy for a run with homogeneous forcing function 
(solid line) and a forcing function whose helicity varies sinusoidally throughout the 
domain (dotted line) simulating the effects of equators at the two nodes of the sinusoidal 
helicity profile. 



A somewhat surprising property of the models with variation of helicity in 
the z-direction is the fact that the mean field varies mostly in the cc-direction, 
i.e. it follows the variation of the background model only weakly; see Fig. |[ 
Therefore, the mean field must be allowed to be two-dimensional, i.e. 

B(x,z,t) = jBdy/Jdy. (20) 

Similar behavior was also found in simulations with boundaries, especially when 
the aspect ratio was large p3| . In the present context we were able to confirm, 
using a two-dimensional mean field dynamo in periodic geometry, that for a oc 
sinfciz the most easily excited mode varies indeed both in x and z\ see Fig. ffl. 
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Fig. 6. Images of the three components of the mean field (averaged over the y- 
direction) for a run with sinusoidally varying helicity in the z-direction. Note that 
the most pronounced component of the mean field is actually B z (x). The large scale 
field is also visible in a y-slice of B z (last panel). 



In the simulations presented so far, boundaries merely tend to reduce the 
final saturation field strength. Thus, the idea to enhance the large scale field by 
small scale losses is not currently supported by simulations. It is quite possible, 
however, that this is simply a consequence of too simple a representation of the 
physical boundary. In the sun, coronal mass ejections are quite vigorous events 
that are known to shed large amounts of helical magnetic fields p8| , p9pc|j3l| | . 
This kind of physics is not at all represented by adopting vacuum or pseudo- 
vacuum (vertical field) boundary conditions, as was done in Ref. J23|. 

We may then conclude that in simulations of large scale dynamos with rela- 
tively simple boundary conditions, open boundaries tend to be more important 
for large scale fields than for small scale fields. Although more realistic bound- 
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Fig. 7. Images of the three components of the mean field for an a 2 dynamo with 
sinusoidally varying Q-effect in the z-direction. The B x and B z components resemble 
those in the direct simulation shown in Fig. ^| 

ary conditions still need to be considered, it is useful to study more carefully 
whether, on observational grounds, a resistively limited dynamo can indeed be 
clearly excluded. 



5 How long is long? — or what the skin depth has to do 
with the solar cycle 

In this section we want to estimate the amount of magnetic helicity that is to be 
expected for a model of the solar dynamo. We also need to know which fraction of 
the magnetic field takes part in the 11-year cycle. Following an approach similar 
to that of Berger , we can bound the rate of change of magnetic helicity in 
terms of the rate of Joule dissipation, Qjouio, and magnetic energy, M. For an 
oscillatory dynamo, all three variables, H, M , and <5j ou ie vary in an oscillatory 
fashion with a cycle frequency w of magnetic energy (corresponding to 11 years 
for the sun - not 22 years), so we estimate |diJ/d£| Si u\H\ and Qjouic Ss wM, 
which leads to the inequality [§3|Q 

\H\/(2imM) < 4kin, (21) 

where 4kin = (2ti/uj) 1 ^ 2 is the skin depth, here associated with the oscillation 
frequency u). Thus, the maximum magnetic helicity that can be generated and 
dissipated during one cycle is characterized by the length scale \H\/ (2/ioM), 
which has to be less than the skin depth -4km- 

For r\ we have to use the Spitzer resistivity which is proportional to T~ 3 / 2 
(T is temperature), so 77 varies between 10 4 cm 2 /s at the base of the convection 
zone to about 10 7 cm 2 / s near the surface layers and decreases again in the solar 
atmosphere. Using u> = 27r/(llyr) = 2 x 10~ 8 s -1 for the relevant frequency at 
which H and M vary we have i'skin ~ 10 km at the bottom of the convection 
zone and i'skin ~ 300 km at the top. 

This needs to be compared with the value \H\/(2fioM) obtained from dy- 
namo models. Although mean-field theory has been around for several decades, 
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the helicity aspect has only recently attracted significant attention. In the pro- 
ceedings of a meeting devoted specifically to this topic p9| , magnetic helicity 
was discussed extensively also in the context of mean-field theory. However, the 
precise amount of magnetic helicity relative to the magnetic energy, and the pos- 
sibility of helicity reversals at some length scale were not addressed at the time, 
although the the evolution of the current helicity has already been investigated 
in the context of a mean field model |5(J . 

For a sphere (or a half-sphere) with open boundary conditions and volume 
V (for example the northern hemisphere), one has to use the gauge-invariant 
relative magnetic helicity of Berger & Field Q , 

H= f (A + A P ).(B-Bp)dF, (22) 
J v 

where Bp = V x Ap is a potential field used as reference field that has on the 
boundaries the same normal component as B. Any additional gradient terms, 
V(/>, in A or Ap yield only a surface term, 

0(B-B P )-dS, (23) 

dV 

which vanishes because Bp • n = B • n. When applied to an axisymmetric mean 
field, which can be written as B = b<p + V x acfi, it turns out that the relative 
magnetic helicity integral is simply |2j| 

H = 2 abdV. (24) 
Jv 

In order to see how the condition ( pl| ) is met by mean-field models, we have 
calculated a typical afl model relevant to the sun and evaluated Eq. ( |24| ) over 
the volume of the northern hemisphere p3[ . 

We recall that in the Babcock-Leighton approach it is mainly the latitudi- 
nal differential rotation that enters. We also note that, although the latitudinal 
migration could be explained by radial differential rotation, meridional circula- 
tion is in principle able to drive meridional migration even when the sense of 
radial differential rotation would otherwise be wrong for driving meridional mi- 



gration [[15 51,52]]. Therefore, we have considered the simple models in spherical 
geometry. The results of such a model show that once B po ] c / B\, c \ t 1S m * ne 
range consistent with observations, B po \ c / Bbcit = (1...3) x 10 -4 , H-^/(2^,qM-^R) 
is around (2 — 5) x 1CP 4 for models with latitudinal shear. (Here, the subscript 
'N' refers to the northern hemisphere.) This confirms the scaling 

H N /{2^M N R) = 0{B pol /B toI ) S> Bpoie/Steif (25) 

Given that R = 700 Mm this means that iJ N /(2/i M N ) « 70. ..200 km, which 
would be comparable with the value of £ s kin near the upper parts of the solar 
convection zone, or for models with only latitudinal shear. 

The surprising conclusion is that the amount of mean field helicity that needs 
to be generated in order to explain the large scale solar magnetic fields is so small, 
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that it may be plausible that microscopic magnetic diffusion could still play a 
role in the solar dynamo. In other words, open boundary effects may well be 
important for understanding the time scale of the dynamo, but the effect does 
not need to be extremely strong. 



6 Connection with a-quenching 

Nonlinear helical dynamos in a periodic domain are particularly simple. They 
provide therefore an ideal benchmark for models of a-quenching. Before applying 
models to more complicated situations, they should pass the test of predicting 
the right behavior in the simple case of an a 2 -dynamo. By "right behavior" we 
mean that the large scale magnetic field saturates as seen in Eq. @. While 
this is already a r elat ively stringent test that allows us to eliminate some earlier 



models (see Sect. EL2), we point out that the complete time evolution (including 
the early kinematic exponential growth phase) can only be described correctly 
by an explicitly time-dependent evolution equation for a that ensures that the 
magnetic helicity equation is obeyed exactly at all times. 



6.1 The lorentzian quenching formula 

We first assume that both a-effect and the turbulent magnetic diffusivity, rjt, 
are being affected by the magnetic field. In the first class of models we assume 

a = a q(B), ry t = ?7to £?(B), (26) 

i.e. we postulate the existence of an algebraic quenching formula for both a and 
rjt- The models that work best are those with a lorentzian quenching formula, 

q(B) = ^ (lorentzian formula), (27) 

1 + aB /B 2 q 

where a is a dimensionless coefficient and B cq is the equipartition field strength 
with BIJuq = (pu 2 ). 

For an a 2 -dynamo in a periodic box, the coefficient a can be readily deter- 
mined. This is because here we actually know the final saturation field strength; 
see Eq. @. In a mean-field model, on the other hand, where quenching is the 
only nonlinearity, the saturation field strength is proportional to aT 1 / 2 . In or- 
der to work out the coefficients, we consider the a 2 mean-field equation for a 
Beltrami field with wavenumber k\. The evolution of the mean-squared field 
strength is then governed by (l^] 

d g2 ) = J^-^__ 
2 & l + a(B 2 )/B 2 q 

In the steady state, the right hand side of Eq. fl28| ) has to vanish, so 

aoh - Vwk 2 



1 + a(B 2 )/BL 



Tjk( = 0, (29) 
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or, using Eq. (|10[), 

a = a o ~ Vtoh = A 
r]k{ rjkfki 

(see Ref. [jl9| ) . The main point here is that the parameter a scales with the 
magnetic Reynolds number, so it is very large under astrophysical circumstances. 
Moreover, if the growth rate A scales with the inverse turnover time on the 
forcing scale, then the k{ in the denominator cancels and a is proportional to the 
magnetic Reynolds number [^)3|j54| , |55|] . This is also consistent with the analysis 
of Ref. || . 

A consequence of the large value of a is a very slow saturation phase and 
not necessarily a very low saturation level, as is usually believed |3. That the 
saturation phase is slow can be seen, for example, by considering the magnitude 
of the right hand side of Eq. ( p8| ) near saturation, so we put 

(B 2 )/i? c 2 q = (l-6)fc f /fc 1 , (31) 
which means that we are a fraction e away from full equipartition. We then have 

rhsofEq. (28) = Q " fcl V kf. (32) 

Since the field is already strong, the '1+' in the denominator may be neglected, 
so we obtain 



rhs of Eq. (28) 



Qpfci - r] W kf 

(aofei - ??tofc?)(l - e) 



T]k\ = r\k\ ks er\k\. (33) 

1 — e 



The rhs of Eq. (J2q) can be regarded as a local growth rate if the field was 
unchanged. The extrapolated saturation time would therefore be (ery/c 2 ) -1 , which 
is more than a resistive time! 

6.2 Other quenching formulae 

By contrast, consider a quenching of the form 

q(B) = 1 - aB 2 /B^ q (quadratic formula). (34) 
To match the right saturation field strength we have to have 

(a ki - r] W kl)(l - ak{/kx) - r\k\ = 0. (35) 

This gives 

a = -=- > — lor small 77. 



aoh - T] to k{ k f k{ 



(36) 



Again, we calculate the instantaneous growth rate for a field that was close to 
final saturation, i.e. we put (B }/B% = (1 — e)k{/ki. This gives 



rhs of Eq. (28) = (ctofci — ?7to&i) 



2 



aoh-riTokf 
1 ; TT\ l - e ) 



- Vkf (37) 
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Fig. 8. Saturation behavior for the quadratic quenching formula (solid line) compared 
with different lorentzian quenching formulae with different values of _R m = \/(r/ki) 
(broken lines). Note the slow saturation behavior, consistent with Eq. (|l4|). The results 
for the quadratic quenching formula are independent of the value of R m , which is 
therefore inconsistent with Eq. (114). 



or, after some simplifications, 

rhs of Eq. (28) = e(a fci - J/rofci). (38) 

This means that it will only take a few dynamical time scales before the extrap- 
olated field will reach final saturation, which is of course incompatible with the 
magnetic helicity constraint (|l4|). Figure |8| shows the great discrepancy between 
these two quenching formulae. The quadratic quenching formula gives, for dif- 
ferent values of R m = X/(i]kf), always the same saturation behavior, while the 
lorentzian formula gives a more prolonged saturation phase as R m is increased. 

One may still be tempted to expect that there could be many other quenching 
formulae that might work as well. One clear counter example is shown in Fig. ^| 
where we compare the results from a cubic quenching formula with the magnetic 
helicity constraint in Eq. (|l4|). For the cubic quenching we just used the very 

simple formula q = 1/(1 + /3 3 ), where /3 2 = a(B )/B^ q . The departure between 
cubic is not very strong, but clearly noticeable. We note that if we replaced 
r\ — > 1 77, the helicity constraint would actually fit, but of course 77 is an input 
parameter, so we cannot just adopt a different value in the analysis. We may 
therefore conclude that cubic quenching can be ruled out. 

Finally, we consider a quenching formula that has the correct strong field 
asymptotics, q — > /3 -2 , and also the same weak field Taylor expansion as the 
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Fig. 9. Saturation behavior for cubic quenching (solid line) compared with the mag- 
netic helicity constraint (dashed line). Note also that if we used the wrong microscopic 
diffusivity, rj — > |tj, the helicity constraint would actually fit. 



lorentzian formula, /3 



1 /„ tan" 1 x /3/9 2 \ , „, ,—2, 
J 2 



q (B) = ± 1 - V ^ j where p 2 = a{ ^)/Bl v (39) 



This expression entered in the quenching formula derived by Field, Blackman, 



& Cho |57| . In Fig. 10 we compare the saturation behavior for the atan formula 
with the lorentzian one. Clearly, the atan formula is much better than any of the 
other formulae considered in this subsection, but the lorentzian formula is still 
considerably closer to the magnetic helicity constraint than the atan formula. 



6.3 Non-universality of the lorentzian quenching formula 

We may now be under the impression that the lorentzian formula is probably the 
correct quenching expression. While is does indeed provide a good description of 
what is going on in the simulations, we note that there are also a few problems. 
Firstly, when applied to other models where the field is in general no longer 
isotropic and force- free (e.g., if there is shear or if there are boundaries) the best 
fit value of a is no longer equal to the value calculated for the a 2 -dynamo; see 
Eq. ©. 

The other problem is rather an uncomfortable prediction from the quench- 
ing model with the lorentzian formula. Since both a and % are quenched by 
equal amounts, % will become comparable with r\ near full saturation. As a con- 
sequence, the cycle period of a.T2-dynamo models becomes comparable to the 
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Fig. 10. Comparison of the saturation behavior for the lorentzian and atan formulae, 
together with the corresponding fits obtained from the helicity constraint. 



resistive time which is rather long. So, mean- field theory of the solar dynamo 
may face a very serious problem, unless open boundary effects (Sect. |]) play an 
important role. Simulations perhaps seem to point into this direction as well: the 



cycle period found in simulations with shear 58 was already rather long, and 
with a further reduction of the magnetic diffusivity by a factor of 2.5, the cy- 
cles disappeared altogether |2j|. On the other hand, the absence of cycles could 
have been for other reasons, for example due to too restricted a geometry (with 
sinusoidal shear flow on a scale only 5 times larger than the energy carrying 
scale of the turbulence). Thus, these two simulations are perhaps not yet fully 
conclusive. 

In the following we point out that there is yet another possibility that is 
at least equally well in agreement with our a 2 -dynamo benchmark result and 
theoretically more theoretically more appealing because it satisfies the magnetic 
helicity equation exactly at all times. 
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Fig. 11. Dependence of a = £ y /Bo y on the residual helicity, obtained by imposing a 
uniform magnetic field Bo = (0, Bo y , 0) and driving the turbulence either through the 
momentum equation (diamonds) or by an extra forcing term in the induction equation 
(asterisks) . 



7 Dynamical quenching 

A quenching formula that we have not yet discussed is the formula for the residual 
O'-effect. 

T T 

a = a K + a M with a K = - - (u ■ u) , a M = +- — (j-b), (40) 

3 3p 

which is due to Pouquet, Frisch, and Leor at p6| , and has frequently been used 
in connection with a-quenching [ p9pc| , |6l| , |62[ . The result that the a-effect is 
proportional to the residual helicity, (u) ■ u) — (j • b) / po has also been confirmed 
numerically (63) by imposing a uniform magnetic field and driving the turbulence 
either through the momentum equation (as is done in the rest of the paper) or 
through a forcing term in the induction equation; see Fig. 11. 

The question now is how to use Eq. (|4(]) in a mean field model, which only 
knows about the mean field, B. If we were to approximate (j • b) in a direct 
manner by B , this would correspond to the quenching formula (|3~4|), which 
is clearly ruled out because it violates the magnetic helicity equation. On the 
other hand, we can explicitly make sure that the magnetic helicity equation (^) 
is obeyed. The contribution (A • B) from the large scale fields to the magnetic 
helicity equation is automatically taken into account by the mean-field equation, 
so we only need to solve for the missing contribution from small scales, (a • b), 
and the two equations for (A • B) and (a • b) must be fully coupled. 

Another way of seeing this |6J] is that, while (A ■ B) stays close to zero 
on short enough time scales, any increase of the magnetic field by the a-effect 
leads to an increase of the large scale magnetic helicity, (A • B). This can only be 
consistent with an almost unchanged (A-B) if there is a simultaneous generation 
of small scale magnetic helicity, (a • b) , of opposite sign, so that 



(A • B) = (A • B) + (a • b) 



(41) 
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stays close to zero. The price to pay for this is that the small scale magnetic 
helicity can also produce an a-effect, om = ^ T (j ' b}/poj but it has the opposite 
sign than ok, so the residual a-effect becomes quenched. Mathematically, this 
quenching of a can be described by the magnetic helicity equation. The contri- 
bution of the large scale field to the magnetic helicity equation follows from the 
mean-field equation. The small scale contribution is exactly such that the sum 
of these two equations gives Eq. (|?]). Thus, we have a pair of two equations Q 

^(A-B)=2(5-B)-2 W o(J-B), (42) 
^(a-b)=-2(£-B)-2 Wo (j-b). (43) 



where £ = u x b is the mean turbulent electromotive force, for which we adopt 
the usual mean-field closure in terms of a-effect and turbulent magnetic diffu- 
sivity, i.e. 

£ = aB 77 tA i J. (44) 

Note that Eq. ( f42|) follows directly from the usual mean-field dynamo equation 
(jjl). Making use of the relation ^ (j • b) = kf (a • b) in Eq. ( [Io| ) , Eq. (|||) becomes 

^ = -2^(i? m ^+a M ), (45) 

where R m = rjt/rj is the appropriate definition in the present context. 

An evolution equation for a was already proposed twenty years ago 
see also Refs |66| , |67|| . Nevertheless, dynamical quenching was usually ignored, 
although it has sometimes been used in mean-field models with the main moti- 
vation to promote and study chaotic behavior in stellar dynamos p^ , |69|j70| j7l|| . 
Kleeorin et al. |2^| were the first to point out that the catastrophic quenching 
of Vainshtein & Cattaneo |54| is just a special case of dynamical quenching. 



7.1 Adiabatic approximation and force- free degeneracy 

Near the saturated state the explicit time derivative in Eq. ([45|), daivi/di, can 
be neglected and the value of cxm adjusts 'adiabatically' as the field saturates. 
Thus, we have ||(J 

= i? m ^-^ + a M , (46) 
or, after substituting ayi — a — ax, 

R m (a(B 2 ) - nt(S ■ B>) + (a - a K ) = 0, (47) 



which yields 

aK + i?m?7t(J • B) 



a 



l + i?m(B 2 ) 



(48) 
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Fig. 12. The early saturation phase of the energy of the mean field in the dynamical 
quenching model for four different values of the magnetic Reynolds number. 



This equation was already obtained by Gruzinov & Diamond [ p9|]67| ] . The late 
saturation phase of a 2 -dynamos is well described by Eq. (Esh . This is because 
near saturation time dependence is governed by the slow resistive adjustment in 
the mean field equation for the large scale field, whilst the a equation is quickly 
adjusting to whatever the large scale field is at any time. 

The reason why the lorentzian quenching formula describes the resistively 
limited quenching behavior so well is because in the case of a nearly force free 
large scale magnetic field the lorentzian quenching formula and the adiabatic 
approximation become identical. Indeed, if the large scale magnetic field is force- 
free, we have (J • B)B = (B")J, which allows us to write the full electromotive 
force, £ = cvB — rjtfioJ , in the form 



a K + R mVt ii (J • B)/B 2 

£ = —o B - ?/toMoJ 

l + R m {B )/B e 2 q 

a K B ?7toMoJ 



- , (49) 

l + i? m (B )/£J2 q l + R m {B)/B^ 

which shows that in the force free case the adiabatic approximation together 
with constant (unquenched) turbulent magnetic diffusivity becomes equal to the 
pair of expressions where both a and rj t are catastrophically quenched. This is 
called the force-free degeneracy ]56| . This degeneracy is lifted in cases with shear 
or when the turbulence is no longer fully helical. 



7.2 a 2 -dynamos 

When applied to an a 2 -dynamo model with dynamical quenching, the helicity 
constraint is well satisfied (5(||64| and the difference between the solutions with 
dynamical and algebraic quenching turns out to be small if R m is less than 
about 1000. The difference in the evolution of magnetic energy with dynamical 
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Fig. 13. The effects of assuming % oc a in the dynamical and the lorentzian quenching 
models for four different values of the magnetic Reynolds number. 
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Fig. 14. Comparison of ai?-dynamo models with algebraic (solid line) and dynamical 
quenching (grey line). f2'/(rjrkf) = 100, aK/Ofrfci) = 0.1, Rm = rjt/v = 10. 



quenching and with algebraic (or lorentzian) quenching, is only a few percent if 
the magnetic Reynolds number is small. The difference increases as the magnetic 
Reynolds number increases; see Fig. [l3[ where we plot the evolution of (B ) 
around the time when the kinematic exponential growth turns into the resistively 
limited saturation phase which was already described in Sect. 3.5. Conclusive 
agreement with simulations is at this point not possible, mostly because the 
magnetic Reynolds numbers are not large enough. 



7.3 aJ?-dynamos 

When shear is included, toroidal field can be regenerated solely by the shear 
term. This is where dynamical and algebraic quenching lead to very different 
behaviors. With algebraic quenching, the reduction of rjt leads to resistively 
long cycle periods in the nonlinear regime. With dynamical quenching, r]t is 
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Fig. 15. Saturation energy versus R m for « 2 -dynamos with boundaries. Models with 
(q — 1) and without (q = 0) loss term are indicated by open circles and open squares, 
respectively. The line gives the fit (B 2 > ~ R~°- m . Simulations of Ref. @ are shown 
as full circles. 



constant and the cycle frequency remains of order unity; see Fig. Here, R m 
is only 10, but for larger values the final field amplitude and the cycle period are 
considerably enhanced. 



7.4 Open boundaries 

Finally we consider a model with vertical field boundary conditions (B x = B y = 
0onz= ±7r) and solve the cv 2 -dynamo equation. The resulting mean-squared 

field strength is plotted in Fig. [inversus R m . We find that (B 2 ) ~ i?" 1 , which is 
consistent with Eq. (|16|), but steeper than what was obtained in the simulations 



| {20| . Kleeorin et al. |22| , [72| pointed out that in Eq. (45) there should be an 
additional loss term on the right hand side. By making this loss term suitably 

2 

i? m -dependent, one could in principle make the i? m -dependence of (B } less 
steep, but some reduction is already obtained by allowing for a diffusion-like loss 
term of the form g?7TV 2 aM on the right hand side; see the open circles in Fig. |l5|. 
The parameter q is used to regulate the efficiency of this loss. The simulations 
of Ref. [EOfl are shown as full circles. 



7.5 Generalization to nonuniform a 

In all astrophysical bodies there are opposite signs of kinetic helicity, (u) ■ u), 
in the northern and southern hemispheres. In the approach of Kleeorin and 
collaborators ]65| , |66| , |67| the dynamical a-quenching framework was always used 
as a theory for nonuniform a, which does not seem to be well justified. Most 
importantly, the connection between (j-b) (in the expression for qm) and (a-b) 
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Fig. 16. Comparison of the horizontally averaged kinetic and current helicity densities, 
(w ■ u) and (j ■ b) , respectively (left hand panel) , and the horizontally averaged magnetic 
helicity densities in the = gauge (a • b) and the Coulomb gauge (a c ■ b) . 



is no longer straightforward when the angular brackets denote ensemble averages 
(which are not really of practical interest) or averages over one or two periodic 
coordinate directions. Indeed, (a-b) is no longer gauge-invariant, so one has to fix 
the gauge. The Coulomb gauge is the most common one, but one should realize 
that going to another gauge can make a major difference. In Fig. |l6| we show 
the magnetic helicity in the cf> — and the Coulomb gauges, (a • b) and (a c • b), 
respectively (see Sect. ||). Note that the Coulomb gauged magnetic helicity is 
about twenty times smaller than that in the <fi = gauge. More importantly, the 
magnetic helicity is not a positive multiple of the current helicity, (j ■ b). (We 
recall that for homogeneous turbulence, (j • b) = kf(a ■ b).) Another problem 
is that when solving numerically the evolution equation for a space-dependent 
a-effect one needs for reasons of numerical stability a diffusion term [7^j7^] . 
However, it is now clear that a loss or exchange of small scale helicity leads to 
an enhancement of the large scale field |5(| , but from simulations we know that 
the presence of an equator rather lowers the energy of the mean field. Perhaps 
this could be fixed by adopting an additional loss term in the mean-field equation 
for the large scale field, but this procedure would be completely ad hoc. One may 
hope however that some kind of a generalization of the dynamical a-quenching 
is at least in principle possible. 



8 Conclusions 

In this review we have outlined some of the main results of isotropic MHD sim- 
ulations in the presence of helicity. We have focussed on the connection with 
the a-effect in mean-field dynamo theory. We should emphasize that in the case 
where the magnetic energy density is uniform in space, the agreement between 
simulations and theory is now well established. In all other cases, things are im- 
mediately more complicated. Moreover, dynamical quenching cannot readily be 
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generalized to the case where qm varies in space. In that case the equation for the 
magnetic helicity density would not be gauge-invariant. Another problem arises 
when qk varies in space and if it changes sign across the equator, for example. 
These are very important aspects requiring clarification. It is quite possible that 
significant improvement in the theory will soon be possible. Without a corre- 
sponding generalization of dynamic a-quenching, if it is ever possible, it would 
be difficult to use dynamo theory in astrophysically interesting circumstances. 
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